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VHE Gamma Rays from PKS 2155-304 
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ABSTRACT 

The close X-ray selected BL Lac PKS 2155-304 has been observed using the Uni- 
versity of Durham Mark 6 very high energy (VHE) gamma ray telescope during 1996 
September/October /November and 1997 October /November. VHE gamma rays with 
energy > 300 GeV were detected from this object with a time-averaged integral flux of 
(4.2 ± 0.7 s tat i 2.0 S?/S ) x 10~ n cnr 2 s _1 . There is evidence for VHE gamma ray emis- 
sion during our observations in 1996 September and 1997 October/November, with the 
strongest emission being detected in 1997 November, when the object was producing 
the largest flux ever recorded in high-energy X-rays and was detected in > 100 MeV 
gamma-rays. The VHE and X-ray fluxes show evidence of a correlation. 



Subject headings: BL Lacertae objects: individual (PKS 2155-304) 
— gamma rays: observations 



galaxies: active 



1. Introduction 

Close X-ray selected BL Lacs (XBLs) are sources of very high energy (VHE) gamma rays at 
energies above 300 GeV. The BL Lac first detected as a source of VHE gamma rays was Mrk 421 
( Punch et al. 1992), following its discovery as a GeV source by the EGRET experiment on the 
Compton Gamma Ray Observatory ( Lin et al. 1992| ). Mrk 421 has been extensively monitored in 
VHE gamma rays and exhibits complex behaviour, having a low-level quiescent state with flaring 
on time scales as short as 30 minutes (Gaidos et al. 1996 ; [Zweerink et al. 1997 ; Takahashi et al. 



1998 ). Mrk 501 is also a source of VHE gamma rays ( Quinn et al. 1996 ), although not detected 



at GeV energies with EGRET, and also exhibits low- level emission with flaring ( Catanese et al. 



1997a| ). In 1997, extremely strong outbursts of TeV emission were detected lasting several months 



(peckers et al. 1997| , IQuinn et al. 1997| , payashida et al. 1998[) The BL Lac 1ES 2344+514 also 
emits VHE gamma rays, most of the evidence for emission coming from a single outburst flCatanese 
et al. 1997b , Catanese et al. 1998 ). All of these objects are close (z ~ 0.03), X-ray selected BL 
Lacs. There have been no reported detections of VHE gamma rays from radio-selected BL Lacs 
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(RBLs — see e.g. Kerrick et al. 1995| , Roberts et al 1998 ) although this category of object is 
frequently detected at GeV energies by the EGRET experiment. 

I I 

Stecker et al. 1996 have interpreted the gamma ray results in the GeV - TeV range and pro- 
pose a model in which RBLs will be GeV gamma ray sources and XBLs will be TeV sources. They 
associate the emission of TeV 7-rays from XBLs with the presence of relativistic electrons with en- 
ergies higher than those in RBLs. They then go on to show that a simple synchrotron self-Compton 
model can explain the differences observed between RBLs and XBLs provided the attenuation of 
the VHE gamma ray flux by pair production with the intergalactic infrared background is taken 
into account ( Stecker et al. 1992 ; Stecker fc de Jager 1997a)) . On the basis of this model, PKS 
2155-304, despite having a redshift of 0.117, is predicted to be a strong TeV gamma-ray source. 

The BL Lac PKS 2155-304 was discovered as an X-ray source during observations made with 
the HE A 0-1 satellite ( [Schwartz et al. 1979] ; Griffiths et al. 1979| ) at a position where the Ariel V 
satellite had previously detected confused emission (|Cooke et al.~197S| ). In many ways, PKS 2155- 
304 may be regarded as the archetypal X-ray selected BL Lac object; like most BL Lac objects it 
is associated with a compact, flat spectrum radio source and has an almost featureless continuum 
which extends from radio to X-ray energies. It is the brightest known BL Lac at UV wavelengths, 
and the object's maximum power is emitted between the UV and the soft X-ray range (Wandel 



fc Urry 1991 ). PKS 2155-304 has a history of rapid, strong broadband variability and has been 
the subject of several multiwavelength monitoring campaigns (see e.g. Brinkmann et al 1994 ; 



Courvoisier et al. 1995| ; Pesce et al. 1997 ). In 1997 November, contemporaneous with some of the 



observations reported here, X-ray emission was detected with the Beppo-SAX satellite flChiappett: 



& Torroni 1997) with a flux equal to the strongest previous outburst. 



The EGRET experiment on board the Compton Gamma Ray Observatory was used to detect 30 
MeV - 10 GeV gamma ray emission from PKS 2155-304 during 1994 November 15 - 29 QVestrand et] 



al. 1995| ). These observations indicated a very hard spectrum, with an integral power-law spectral 
index of 1.71 ± 0.24, and this, combined with its proximity, make it an excellent candidate TeV 
source. However, a recent detailed analysis of the shape of the X-ray spectrum suggests that PKS 
2155-304 may not be an important TeV emitter, producing photons with a maximum energy of 
~ 800 GeV, although a very high degree of beaming may allow the production of photons with 



higher energies ( Giommi et al. 1998 ). Indeed, the absence of any signature of 7-ray absorption due 
to pair production at energies up to 8 GeV in observations of PKS 2155-304 is consistent with the 



suggestion that it contains a relativistic jet oriented nearly along our line of sight ( Vestrand et al. 
19951) . 



2. Observations 
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2.1. The Mark 6 Telescope 

The University of Durham Mark 6 atmospheric Cerenkov telescope has been in operation at 
Narrabri, NSW, Australia since July 1995. The telescope is described in detail elsewhere ( [Arm- 
strong et al. 199~7| ). It uses the well-established imaging technique to separate VHE gamma rays 
from the cosmic ray background, combined with a robust noise-free trigger based on the signals from 
three parabolic flux detectors of 7 m diameter and aperture f/1.0 mounted on a single alt-azimuth 
platform. A 109-element imaging camera (0.25° pixel size) is mounted at the focus of the central 
flux detector, with low resolution cameras each consisting of 19 pixels (0.5° pixel size) mounted at 
the focus of the outer (left and right) flux collectors. The telescope is triggered by demanding a 
simultaneous temporal (10 ns gate) and spatial (0.5° aperture) coincidence of the Cerenkov light 
detected in the three cameras. This triggering system enables the telescope to detect low energy 
gamma rays with high immunity from triggering by local muons. Initial work suggests that our 
system has a detection probability for gamma rays of about 1% at 100 GeV with the probability 
of detection rising slowly to about 40% at 250 GeV. 

The Mark 6 telescope has been designed to provide stable operation which allows the obser- 
vation of weak DC sources. The multiple-dish triggering system is stable against variations in 
performance due to accidental coincidences. All detector packages are thermally stabilised. The 
atmospheric clarity is continuously monitored both using a far infra-red radiometer and an axial 
optical CCD camera which enables the position and visual magnitude of guide stars to be moni- 
tored. The gain and noise performance of the PMTs, digitizer pedestals and associated electronics 
are continuously monitored by: 

1. triggering the telescope at random times using a nitrogen laser / plastic scintillator / optical 
fibre light guide / opal diffuser system to simulate Cerenkov flashes and so enable flat-fielding, 
and 

2. producing false triggers at random times to measure samples of the background noise. 

PKS 2155-304 was observed in 1996 September/October /November and 1997 October/November 
under moonless, clear skies. The 1997 November observations were made as part of a multiwave- 
length campaign co-ordinated by T. Vestrand. An observing log is shown in Table |]. Data were 
taken in 15-minute segments. Off-source observations were taken by alternately observing regions 
of sky which differ by ± 15 minutes in right ascension from the position of PKS 2155-304 to en- 
sure that the on and off segments possess identical azimuth and zenith profiles. This off-source 
- on-source - on-source - off-source observing pattern is routinely used to eliminate any first or- 
der changes in count rate due to any residual secular changes in atmospheric clarity, temperature 
etc. The gross counting rate of the telescope was typically ~ 800 per minute near culmination 
(zenith angle 5°) reducing to ~ 150 per minute at zenith angles of 60°. Most of the observations 
at large zenith angles were made in 1997 November, in order to maximize exposure during the 
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Date 




No. of 


Date 




No. of 






scans 






scans 






ON source 






ON source 


1996 


September 3 


2 


1997 


October 20 


3 


1996 


September 5 


7 


1997 


October 22 


7 


1996 


September 6 


5 


1997 


October 23 


7 


1996 


September 7 


6 


1997 


October 25 


3 


1996 


September 8 


4 


1997 


October 29 


4 


1996 


September 10 


8 


1997 


October 30 


2 


1996 


September 12 


3 


1997 


November 18 


2 


1996 


September 14 


7 


1997 


November 19 


4 


1996 


September 15 


7 


1997 


November 20 


5 


1996 


September 17 


5 


1997 


November 22 


2 


1996 


September 18 


3 


1997 


November 23 


4 


1996 


October 9 


4 


1997 


November 24 


4 


1996 


October 10 


8 


1997 


November 25 


2 


1996 


October 12 


5 








1996 


October 13 


4 








1996 


October 14 


6 








1996 


November 2 


3 








1996 


November 7 


6 








1996 


November 10 


4 








1996 


November 11 


5 








1996 


November 12 


5 









Table 1: Observing log for our observations of PKS 2155-304 during 1996 and 1997. 

multiwavelength campaign. Detections of Mrk 501 during its 1997 flaring activity using the Mark 
6 telescope at zenith angles of > 70° have demonstrated the usefulness of TeV observations using 
the atmospheric Cerenkov technique at very large zenith angles flChadwick et al. 1998c ). 

Data were accepted for analysis only if: 

1. the sky was clear and stable, and 

2. the gross counting rates in each on-off pair were consistent at the 2.5 a level. 

When these criteria are met, we have a total of 41 hours of on-source observations with an 
equal quantity of off-source data. 
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3. Results 

Routine reduction and analysis of accepted data comprises the following steps: 

1. calibration of the gains and pedestals of all 147 PMTs and digitizer electronics within a 15 
minute segment, using the embedded laser and false coincidence events, 

2. software padding of the data to equalize the effects of on- and off-source photomultiplier noise 
on data selection flCawley et al. 1993 ; Fegan 1997 ), 



3. identification of the location of the source in the camera's field of view for each event, using 
the axial CCD camera, 

4. a calculation of the spatial moments of each shower image relative to the source position, and 

5. rejection of events containing an image which would be unlikely to be produced by gamma 
rays. 

Events considered suitable for analysis are those which are confined within the sensitive area 
of the camera (i.e. within 1.1° of the centre of the camera) and which contain sufficient information 
for reliable image analysis, i.e. which have SIZE > 500 digital counts, where 3 digital counts ~ 1 
photoelectron, and 200 digital counts are produced by a 125 GeV gamma ray. 

Monte Carlo simulations indicate that the image shape of the Cerenkov light from a 7-ray 
shower can be approximated by an ellipse, the major axis of which is oriented towards the source 
position, whereas a cosmic ray (hadronic) shower produces a broader, more irregularly shaped 
image flHillas 1985| ). The image can be parameterized using techniques developed by the Whipple 



group which describe both the shape and the orientation of the image. In addition a measure of the 
fluctuations between the centroids of the samples recorded by the left and right flux collectors of 
the Mark 6 telescope provides a further discriminant, -Ddist ( phadwick et al. 1998a ). Gamma rays 



are identified on the basis of image shape and left/right fluctuation, and then plotting the number 
of events as a function of the pointing parameter ALPHA; 7-ray events from a point source will 
appear as an excess of events at small values of ALPHA. 

The selection criteria applied to these data are summarized in Table 0. They constitute a 
standard set of criteria developed to include allowance for the variation of parameters with event 
size and are routinely applied to data from all objects recorded at zenith angles less than 45°. Small 
changes in the values for WIDTH, CONCENTRATION and £> d ist are made for data recorded at 
zenith angles greater than 45°. These selection criteria differ little from those used in our earlier 
reports of the detection of VHE gamma-rays from PSR B1706-44 QChadwick et al. 1998b|) and 



Cen X-3 (Chadwick et al. 1998a); a decrease in the weight given to the signal recorded from the 



outer (guard ring) detectors in the camera has produced an A LPH A-plot for the background with 
an increasing population of large ALPHA events. 
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Parameter Ranges Ranges Ranges Ranges Ranges 

SIZE (d.c.) 500 - 800 800 - 1200 1200 - 1500 1500 - 2000 2000 - 10000 

DISTANCE 0.35° - 0.85° 0.35° - 0.85° 0.35° - 0.85° 0.35° - 0.85° 0.35° - 0.85° 

ECCENTRICITY 0.35 -0.85 0.35 - 0.85 0.35 -0.85 0.35 -0.85 0.35 - 0.85 

WIDTH < 0.10° < 0.14° < 0.19° < 0.32° < 0.32° 

CONCENTRATION < 0.80 < 0.70 < 0.70 < 0.35 < 0.25 

D dist < 0.18° < 0.18° < 0.12° < 0.12° < 0.10° 



Table 2: The image parameter selections applied to the PKS 2155-304 data recorded at zenith 
angles less than 45°. For data recorded at zenith angles greater than 45° small modifications are 
made to the WIDTH, CONCENTRATION and D dist criteria. 

The number of events remaining ON and OFF source after the application of the selections 
described above are summarized in Table ^. The ALPHA distribution for the whole dataset, after 
subtracting the results for OFF from ON observations, is shown in Figure [l], together with the 
ALPHA distributions for the ON and OFF source data individually. No normalization of ON and 
OFF data rates has been applied. There is an excess of events at small ALPHA, the expected 7-ray 
domain, and imposing a selection of ALPHA < 22.5° yields a gamma ray detection significant at 
the 6.8 a level for the total dataset. 



On Off Difference Significance 



Number of events 1021083 1023280 -2197 -1.5 a 

Number of size and 600856 598733 2123 1.9 a 
distance selected events 

Number of shape 37125 36151 974 3.6 a 
selected events 

Number of shape and 6099 5370 729 6.8 a 
ALPHA selected events 



Table 3: The results of various event selections for the PKS 2155-304 data. Data from observations 
at all zenith angles have been combined. 

In addition, a false source analysis has been performed by re-analysing the ON and OFF- 
source data using a matrix of assumed source positions in celestial co-ordinates. For each trial 
source position in the matrix, DISTANCE and ALPHA are recalculated for each event, and a 
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Fig. 1— (a) The ALPHA distributions ON and OFF source for PKS 2155-304. The grey line 
refers to ON source data, (b) The difference in the ALPHA distributions for ON and OFF source 
events. 



selection of events is made with values of DISTANCE as above and ALPHA < 22.5°. The result is 
shown in Figure ^; the excess of gamma ray candidates maximises when the source is assumed to 
be at the centre of the plot. The central point in the matrix is the position of PKS 2155-304; it is 
not the geometrical centre of the camera, which is typically 0.1° — 0.2° from the source position. 



3.1. Observed flux 



544 excess events identified as 7-rays were detected in 32.5 hours of on-source observation at 
zenith angles less than 45°. The current selection procedure for data recorded at zenith angles 
< 45° is estimated to retain ~ 20% of the original gamma ray events. The collecting area has 
been estimated from Monte Carlo simulations and is ~ 5.5 x 10 8 cm 2 for these observations. Using 
this estimate, the 7-ray flux incident on the earth's atmosphere was (4.2 ± 0.75 s t a t i 2.0 sys ) x 
10~ n cm~ 2 s _1 . The average gamma ray threshold for this observation, defined as the energy at 
which the trigger probability is ~ 50%, was 300 GeV. 

We also note that the high energy (large size) bin for our observations at zenith angles greater 
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Distance from source (deg) 



Fig. 2. — The significance of the number of excess events as a function of assumed source position 
for PKS 2155-304. The grey scale is such that black corresponds to a detection probability of 
> 6 a. Contours are at 0.7 a intervals. 

than 45° yields a detection significant at the 2a level. We estimate that this sample corresponds 
to a minimum 7-ray energy of about 3 TeV and so there is evidence for detection of photons with 
energies > 3 TeV from PKS 2155-304. 



3.2. Time variability 

A characteristic of previous observations of VHE gamma rays from XBLs has been the strong 
time variability. We have investigated our data for variations in the TeV flux from PKS 2155-304 
on timescales of months and days. 

We show in Figure || the variation in the detected 7-ray flux averaged over the observations 
made in the observing periods in 1996 September, October and November and in 1997 October and 
November. The strength of the 7-ray emission is defined as the number of 7-ray candidates (shape 
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Fig. 3. — The variation of the VHE gamma ray flux from PKS 2155-304 averaged over the observing 
periods in 1996 September, October and November and in 1997 October and November. 



and orientation selected events) to the number of background protons during the same observation. 
This makes some allowance for variations in the sensitivity of our telescope as the zenith angle 



changes. Using the test for constancy of emission used by the EGRET group (McLaughlan et al. 



1996), the data suggest that the VHE 7-ray emission is time variable. The largest excursion from 
the average flux occurred during observations in 1997 November. 

The variation of the gamma-ray emission during observations made in individual nights is 
shown in Figure ||. We have no evidence for strong gamma-ray flaring on a time scale of days during 
these observations. However, the data suggest that the VHE gamma-ray emission is variable on a 



timescale of days, using the test of McLaughlan et al. 1996 . 



4. Discussion 

PKS 2155-304 is the fourth X-ray emitting BL Lac to be established as a VHE emitter and, 
thus far, the furthest from earth with a redshift of 0.117. The emission shows the features of time 
variability demonstrated in the other BL Lacs detected at VHE energies. 

In Figure |5| we show the spectral energy distribution (SED) from PKS 2155-304 from radio 
to VHE gamma ray energies, including the present results. The present VHE result is consistent 
with the SEDs for the other VHE-emitting blazars. The VHE behaviour is in general agreement 
with the predictions of Ghiscllini et al. 1998| and, as for the other three established VHE blazars, 



indicates the success of the unified model of AGNs in explaining the gross features of VHE emission. 
However, the present VHE data are unable to distinguish between the predictions for the external 
Compton or the synchrotron self-Compton model, as is the case for TeV observations of the other 
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Fig. 4. — The variation of VHE gamma ray flux from PKS 2155-304 from each observation made 
in (a) 1996 September, October and November and (b) 1997 October and November. 

VHE blazars. 

During our observations of PKS 2155-304 measurements of its 2-10 keV X-ray emission were 
available from the ASM on RXTE. [] In Figure |(a) we show the correlation between the average 
VHE emission during each of the five months when data were available and the corresponding 
average X-ray rate from the RXTE quick look analysis. The correlation coefficient for the data is 
0.66. A similar plot for the daily averaged TeV and X-ray measurements are shown in Figure ||(b). 
In this case, the correlation coefficient is 0.22. 

The data in Figure || show evidence for a correlation between the strengths of VHE gamma 
ray emission and X-ray emission from PKS 2155-304 on timescales of both days and months. In 
particular, the observing periods during which we failed to detect significant VHE gamma-ray 
emission (1996 October and November) coincide with a period of low X-ray emission according to 



the RXTE ASM measurements. We also note that Giommi et al. 1998 based their conclusions that 
PKS 2155-304 should not emit VHE gamma rays from an X-ray observation taken with Beppo-SAX 
shortly after our 1996 November (null) observation. The general correlation between X-ray and 
VHE emission in PKS 2155-304 is in agreement with multiwavelength observations of Mrk 421 and 
Mrk 501 (reviewed by e.g. |Ulrich, Maraschi fc Urry 1997| ). These show such correlated behaviour, 
indicating that the same population of electrons produces the X-ray and TeV emission; this is the 
basis for many models for VHE production in AGNs. 



As has been pointed out by |Catanese et al. 1997a and Dermer 1998] , a detailed study of the 



temporal correlation between the X-rays and TeV gamma-rays during flares can provide a test for 
beaming. However, none of the flaring activity observed in our observations of PKS 2155-304 is at 



Available on the web at http://space.mit.edu/XTE/asmlc/pks2155-304.html. 
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Fig. 5. — The spectral energy distribution of PKS 2155-304. The VHE point is from the present 
work. 

a strength where meaningful limits can be placed on the contribution of beaming. 

The strongest VHE emission from PKS 2155-304 occurred in 1997 November (during the 
multiwavelength campaign). We note that the strongest X-ray emission ever observed from this 
object was detected by Beppo-SAX during observations in 1997 November 22 - 24 ( |Chiappetti 



fc Torroni 1997j ) . EGRET also detected GeV gamma-rays, again at a flux considerably higher 



than prevous detections, during an observation from 1997 November 11 - 17, just before the onset 
of our November observation ( Srcckumar fc Vestrand 1997] ) . The results of the 1997 November 



multiwavelength campaign will be reported elsewhere. 

Our observations indicate that PKS 2155-304 emits VHE gamma rays at energies above about 
300 GeV and that the energy spectrum extends, with significance at the 2a level, beyond 3 TeV. 
The estimated average integral flux above a threshold of 300 GeV is (4.2 ± 0.7 s t a t ± 2.0 SJ/S ) x 
10~ n cnr 2 s -1 . 

Observations of VHE gamma-rays from AGNs can be used to probe the inter-galactic infra-red 
radiation field ( [Stecker et al. 1992 ). The most recent version of this model ( [Stecker fc de Jager 



1997b , Stecker fc de Jager 1998 ) uses estimates of the infra-red radiation field due to Malkan fc 



Stecker 1988 to predict the optical depth for TeV sources as a function of distance. They predict 
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Fig. 6. — The correlation between VHE emission from PKS 2155-304 and the X-ray flux measured 
by the RXTE satellite. Fig 6(a) shows the correlation when averaged over an observing period 
(~ 10 days); the correlation coefficient is 0.66. Fig 6(b) shows the day-to-day correlations; the 
correlation coefficient is 0.22. 



an optical depth of 1 - 2 for an object emitting 3 TeV gamma-rays at z = 0.12. Our detection of 
> 3 TeV gamma-rays from PKS 2155-304 is thus not in conflict with the predictions of Stecker fc 
de Jager 1998 . However, a more sensitive test must await a determination of the VHE gamma-ray 
spectrum from our data. 
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